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SUMMARY 

A detailed method is presented for determining the temperature 
and flow of heated gas necessary for ice prevention of hollow pro- 
peller blades in flight and- icing conditions. The propeller blade 
is analytically divided into a number of short radial segments, 
which are successively treated as separate heat exchangers. Expres- 
sions for the total external and internal heat transfer are com- 
bined to determine the surface temperatures of each segment. The 
thermodynamic steady -flow equation is given for the internal gas- 
flow process and expressions are obtained for the radial variations 
of gas temperature and pressure within the blade. For a given 
initial gas temperature in the blade shank cavity, the minimum gas 
flow is determined, which will provide surface temperatures' of at 
least 32° F everywhere on the heated portion of the blade. 

An expression for the required heat-source input to the gas 
is included and a formula is given for calculating the required 
blade -tip nozzle area. 

A discussion is included of the indicated benefits to be 
derived from certain alterations of the blade internal flow 
passage. 


INTRODUCTION 

Thermal prevention of the formation of ice on aircraft pro- 
pellers has been previously investigated by means of electrical 
and hot -gas heating. The heat requirements associated with the 
use of external, electrically heated rubber blade shoes are 
analyzed in reference 1. Flight investigations of thermal sys- 
tems using hot gas within hollow propeller blades have been 
reported in reference 2, and by Palmatier and Brigham of the 
Curtiss -Wright Corporation. 
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The present analysis was made at the NACA Cleveland labora- 
tory to provide a method for predicting the hot-gas flow and the 
initial gas temperature required for satisfactory ice prevention 
on hollow propeller blades. This method involves a system of 
progressive approximations considering successive radially dis- 
posed blade segments as separate heat exchangers. Account is 
taken of the heat required to raise the temperature of the inter- 
cepted free water to the blade-surface temperature, the cooling 
effect of the evaporation of water, and the kinetic heating of the 
external -air boundary layer.' The physical changes of the hot gas 
in flowing through the hollow blade are determined in order to 
find the required tip nozzle area-, The detailed analysis is 
applied to a typical propeller blade; for flight at two assumed 
operating conditions at the same icing rconditicms to illustrate 
the step-by-step procedure and to demonstrate typical results. 

A method of modifying the blade internal, passage is suggested 
whereby the heating requirements may be: reduced for a given 
application of the hot-gas method : of preventing ice on typical 
hollow propellers. 1 ■ 

. > 

SYMBOLS' • • 

The. following symbols are used in equations taken from 
references : 


A heat-transfer area, (sq ft) 

Ajj cross ^sectional area of blade tip nozzle, (sq ft) 

A p cross-sectional area of blade internal -flow passage, 

(sq ft) 

Cj blade-section lift coefficient 

c blade chord, (ft) 

c_ specific heat of gas for polytropic process, 

” (Btu/(lt>)(°F)) 

c specific heat of gas at constant pressure, (Btu/(lb)(°F)) 

Ir 

D c diameter of cylinder whose radius is equal to leading-. 

c edge radius of given blade section, (ft) 

hydraulic diameter, (4 Ap/P), (ft) 
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e 1-2 

e 

F l-2 

f 

g 

H 

h 

J 

k 

L 

l 

M 

M ev 

m 

Nu 

n 

P 

Pr 

p 

Py 

% 

Qn 

Or 


change of radial kinetic energy per pound of gas in flow- 
ing through radial blade segment, (ft -lb/ lb ) 

base of natural logarithms 

flow energy dissipated by friction' per pound of gas in 
flowing through radial blade segment, (ft -lb/lb) 

friction coefficient , 

acceleration of. gravity, 32.2 (ft/sec 2 ) : . 

rate of heat transfer per unit area, (Btu/(hr)(sq ft)) 

convective- heat -transfer coefficient, (Btu/ (hr)(sq ft)(°F)) 

mechanical equivalent of heat, (ft-lb/Btu) 

thermal conductivity, (Btu/(hr) (sq ft) (°F/ft ) ) 

latent heat of evaporation of water, (Btu/lb), 

radial length of blade segment, (ft) 

rate of interception of water, (lb/(hr)(sq ft)) 

rate of evaporation of water, (lb/(hr)(sq ft)) 

liquid-water content of ambient air, (grams /cu m) 

Nusselt number, (h&j/k or h5 t /k) 

exponent in polytropic process, pv n = constant 

perimeter of blade internal -flow passage, (ft) 

Prandtl number, (3600g) pc^/k 

absolute static pressure, (lb/sq ft) 

pressure of saturated water vapor, (lb/sq ft) 

heat input to gas flow per blade from heat source, 

(Btu/hr) 

heat escape per blade at tip nozzle, (Btu/hr) . 

total heat added per blade to internal gas flo^^rom heat 
source and propeller work, (Btu/hr) 
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^trans 

Q l-2 


<1 

E 


Re D,h . 

Re &,z . 

r 

Srji 


T 

t 

At 

U 

”r 


u t 


V 




Vr 


heat transferred through blade metal from gas flowing 
through radial blade segment, (Btu/hr) 

net heat loss per pound of gas in flowing through radial 

blade segment, Qtrans _ R l-2 ; (Btu/lb) 
w j 

dynamic pressure of ambient air relative to blade sta- 
tion, (lb/sq ft) 

gas constant (for air = 53.3), (ft-lb/(lb)(°E) ) 

Reynolds number based on blade-section chord 

Reynolds number based . on 1 hydraulic diameter 

Reynolds number based on laminar boUridary-layer thickness 

radius at any propeller rblade station, (ft) 

surface total length of heated blade section in chordwise 
direction over both 'camber and -thrust f.ace, (ft) . 

surface distance from blade -section stagnation point to 
any .point on heated surface in chordwise direction, (ft) 

absolute total temperature, (°R) 

static temperature, (°F) 

temperature rise at surface due to boundary -layer 
friction, (°F) 

internal energy of gas, (Btu/lb) 

radial velocity of internal gas relative to propeller 
blade, (ft/sec) - 

tangential velocity of propeller blade at given radius, 
(ft/sec) 

air velocity, (ft/sec) 

average air velocity over camber or thrust face of blade, 
(ft/see) 

resultant or helical velocity of propeller blade at any 
radius, (ft/sec) 
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v specific volume of- internal gas, (cu ft /it) 

W prop increment of enthalpy due to work done by rotating pro- 
peller on pound of gas in flowing through radial blade 
segment, (ft-lb/lb) 

Wi_2 net compression work done on pound of gas in flowing 

through radial blade segment, Wp rop - F^, (ft-lb/lb) 

w rate of internal gas flow per blade, (lb /hr) 

X Hardy's evaporation factor 

x distance from leading edge measured along chord, (ft) 

Z exponent of Pr in determining kinetic temperature rise, 

l/2 for laminar flow and l/3 for turbulent flow 

a angle of attack of blade section, (deg) 

t ' 

7 ratio of specific heats 

6 1 thickness of laminar boundary layer on blade surface, (ft) 

6 t heat -transfer length of turbulent boundary layer, (ft) 

e momentum thickness of boundary layer on blade surface 

at transition, (ft) 

£ turbulent boundary- layer parameter 

turbulent boundary-layer parameter at point of transition 

A propeller drive-shaft torque increment due to internal 

gas flow, (lb-ft/(lb/sec)) 

p absolute viscosity, ( (lb) (sec)/(sq ft) ) 

l) kinematic viscosity,- (sq ft/sec) 

p density, ( (lb) (sec 2 ) /ft 4 ) 

T static temperature, (°R) 

t arithmetic average of static temperatures of ambient air 

.and blade surface, (°E) 

jp angular position relative to stagnation point on blade- 

section leading-edge circular arc, (deg) 
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angle of impingement! .of water droplet on blade surface, 
(deg) 


0) • propeller rotational speed, (radians/sec) 

Subscripts: 

0 ambient atmospheric conditions 

1, 2 internal gas at entering and leaving ends, respectively, 
of a given radial segment , 


av 

b 

d 


g 

i 

m 

s 

w 


external air side 


of . propeller blade 


average 

1 * } 

outer edge of external-air boundary layer 

datum temperature (for determining heat -transfer differ- 
entials) 

final gas conditions (exit from last radial blade segment 
prior, to tip nozzle) 

internal gas side of propeller blade 

initial gas condition (inlet to first radial blade 
segment) 

mean value between points 1 and 2 

external blade surface 

occurrence of condensation or evaporation of water 


Primed symbols denote conditions that are changed because of 
altered blade interiors . 


METHOD FOE DETERMINING DESIGN REQUIREMENTS 
Description of Solution 

The determination of the required gas -flow rates for a hot- 
gas type of anti-icing system on a given aircraft propeller 
involves . a tedious series of calculations. The information 
available at the present time is insufficient to determine * 
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quickly the maximum requirements for a given propeller operating 
through a range of flight and icing conditions. Individual calcu- 
lations must be made for a complete series of selected, critical 
flight and icing conditions. The blade -tip nozzle area and the 
heat-source capacity are then calculated for the maximum require- 
ments so determined. 

The degree. of approximation obtained through the use of the 
subsequent analysis will remain uncertain until more complete and 
reliable experimental data are available on propeller-blade 
external and internal heat -transfer coefficients, kinetic heating 
of wet air, heat losses due to evaporation of water, the extent 
of heat conduction in the blade metal, and the determination of 
minimum standards for ice prevention. 

A schematic diagram of the heatedrgas flow through a typical 
hollow propeller blade is illustrated in figure 1. The gas is' 
first heated by a suitable heat source, enters a stationary 
transfer manifold, and then passes through a collector ring and 
cuffs fastened to the propeller hub and into the blade shanks. 

The hot gas then flows radially outward through the hollow blades, 
which may have internal radial partitions (fig. 1), to discharge 
nozzles in the blade tips. 

The method suggested starts with the selection of a number 
of operating conditions (propeller rotational speed, airspeed, 
density altitude, and ambient-air temperature) characteristic of 
taxiing, take-off , climb, cruise, and maximum speed with the 
determination of the corresponding liquid -water contents from 
the recommendations shown in figure 2. . The ■ required internal 
hot-gas flow is then calculated for each condition using an 
initial _hot -gas temperature as high as structural. considerations 
will permit. 7 ; — — — — 

For some flight conditions the required internal gas flow 
obtained in the calculations using the maximuja initial gas tem- 
perature will be so low that a reduction in the initial tempera- 
ture and an increase in the gas flow might seem advisable. In 
many cases, however, the maximum heat requirements will probably 
be met only by utilizing the highest possible initial hot -gas 
temperature because of limitations of the pr.opeller in pumping 
the gas. The use of high gas temperatures and low gas flows has 
the additional advantage of keeping the pumping power at a 
minimum. - • • 

The required gas flow, which, is. defined, as the rate that will 
maintain the heated portion of a. blade surface at temperatures at 
or above 32° F, must be computed by trial. For conventional 
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propeller-blade designs, the point .of. lowest calculated- surface 
temperature will usually be at the stagnation point of a radial 
station near the hub. Under severe ; icing conditions, there is no 
great variation of calculated leading-edge surface temperature 
with propeller radius at the inboard stations, although heat con- 
ducted radially outward from the. blade hub will actually raise the 
surface temperatures at the inner radii. If the effect of this 
heat conduction from tils blade shank does not extend too far along 
the blade at the leading edge, the use of the stagnation -point 
surface temperature at the . innermost blade segment as the critical 
point for determining the required gas flow will permit more rapid 
progression with the trial- values and Seldom will a complete 
analysis of all segments be necessary';; : , 

' ■ t , 

Steps in Solution 1 • 

For each set of flight and 1 icing conditions the following 
step-by-step method of solution is suggested. 

I. Divide the blade 1 irito^a convenient number of radial seg- 
ments each 2 long. (In the numerical example given 
subsequently, radial segments 1 ft long were employed; 
experience may indicate a need for shorter segments of 
perhaps 6 in. in length.) Tabulate the cross-sectional 
flow areas Ap, the external and internal heat -transfer 

areas A a (= s T l) and A g (= PZ), respectively, the 
diameters of equivalent leading-edge cylinders D c at 
the centers of the segments, the internal perimeters P, 
and the chord lengths c. The cross-sectional flow 
areas should be tabulated for the inlet,, the center, and 
the outlet of each segment. 


II. Determine the External factors that are unaffected by 
the internal heat flow as follows : For the assumed 

propeller-blade speed, find the lift coefficient Cj 
and angle of attack a at the center of each radial 
segment by the method of reference 3 or other accepted 
procedure. With these data calculate the chordwise 
distribution of local air velocity Vj>. 


A. Determine the local .external air-film heat -transfer 
coefficient h a for the laminar regime (from refer- 
ence 4) from 


^a ~ 


Nu k 


( 1 ) 
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where the value. of the laminar boundary-layer thick- 
ness Sj is found from reference 5 as. 

" l s/c 8.17 

m <*> 

0 

For the turbulent regime the heat -transfer coeffi- 
cient has been given in reference 6 as 



•h a = Nhk 
_ a 8. 


(3) 


The value of the turbulent boundary -layer heat- 
transfer length is slightly involved. The 

turbulent boundary^ lay ev. parameter £, as derived 
in reference 7, is related to as follows: 


&t = 


t Z o . 


Re, 


(4) 


% 


The value of the parameter t, denoted at the 

transition point (appendix A), is given in refer- 
ence 7 as 


= 2.56 log e 4.075 


Vb « 


where 


c = 0.289 


(5) 


( 6 ) 


the value of 6^ being taken at the transition 

point. Succeeding values of £' may be found from 
the relation 


d i 


6.15 AV b 
V b 


V v 


fit).. 


dx Vv dx u 
where the function f (£ ) is 

f(£) = 10.411 f 2 g-°.3914t 


(7) 


( 8 ) 
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An alternate and simpler method than that above for 
determining h a (explained in appendix A) is given 

in reference Q as 



for the laminar regime, and 


( 10 ) 


h a - 0-524 
for the turbulent regime, where . 

•*7 • :'i ‘y 

ff _- ' . \ . >, , _ 

V av ~ V R f 1;± t~~~ — ) (ID 

av a \p f 4 cos a J 

\V v "' ’ *r- j 

The + is used f.or the camber surface and the -/ is 
used for the thrust J face. Inasmuch as the values of 
h a in equations (1), (3), (9), and (10) Gannot be 

determined at the stagnation point, determine values 
of • h a near the leading edge (appendix A) from the 
equation of reference 8 


h a = 0.194 (T y ) 


0.49 / V R p, 

[— 


P0 8 
'c 


0.50 


£| 3 ) 

90 / 


( 12 ) 


B. Determine the rate of water interception M (appen- 
dix B) from 


M= ® V R sin (13) 

C.- Determine the local .heat -transfer datum et.r tempera- 
ture t a £ (appendix B) from 

t a,d t b,.w + ^ t a,v ( 14 ) 

for saturated air flowing over a wet surface. If 
the air flowing over, the surface is unsaturated, the 
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local air temperature at^ the outer edge of the 
houndary layer becomes t b , and if the surface is 
dry the air-temperature rise in the boundary layer 
due to friction becomes At a . The values of t b w 

and t b are determined from the adiabatic compres- 
sion lines shown in figure 3, plotted from data of 
reference 9. From the intersection point of the 
ambient -air temperature t 0 and pressure p Q , 

determine the proper, adiabatic line; then the local 
temperature at the outer edge of the boundary layer 
t b at any point- over the airfoil is found by fol- 
lowing this adiabatic line parallel to those shown 
to the 3,dc al pressure p^ given by 




(15) 


The values of At a and At a w have been developed 
in reference 10 as follows: 


and 


At a = 


V Pr 2 

2gJCp 


At 


a,w 


= At n - 0.622 


(- ^ r^) 


(16) 


(17) 


where the saturated vapor pressures 


p v,a 


and 


correspond to the temperatures - t & ^ and t b , 
respectively, necessitating a solution of At a w 
trial at each point, 3 


p v,b 

ty 


III. In order to determine the changes in the internal gas 

flow through the blade, select the desired inlet -gas tem- 
perature tg^i, which equals t g for the first radial 

segment, and a trial value of hot-gas flow w. 

\ 

A. Assume the temperature cf the gas leaving the first 
segment t g g • 
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1. Determine the internal gas -film heat -transfer 
coefficient hg (appendix C) from reference 11 


h g = 4.1 x 10 
using 


(t ) 0,3 w 0,8 P°’ 2 
•4 8} m 


(18) 


T . - -gjJl! ii i 

.8i» • 2 


(19) 


2. Determine the mean heat-transfer datum gas tem- 
perature for the segment tg <j_ from 


1 <u r>m ) 2 Pr°> 33 

t g,d = t g,m + At g te :fc S>m + ’ 2gJc_ " 


(20) 


where 


u 


r,m * ^ 3600 Pg;ia 


_ ( T ) 

w R ' g,nr 


( 21 ) 


For the first estimate of tg ^ 2 , let 


P « p_ i ; for subsequent trials let 

g,m J- 




~ P g>l * P g> 2 


of the preceding trial. (Details for calculat- 
ing Pg 2 presented subsequently.) As 

the solution is approached, this error will 
become very small. 


3. Determine the chordwise distribution of surface 
temperature t s (appendix C) at the center of 


the segment from 


h_t 


*s = 


g^g,* o 


Art + 


k Q Xt Q a + M \tQ + 




g 


a A a,d 
+ .(h a x + M) A a 


R 

2gJ. 


Aft 


( 22 ) 
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whore (see appendix B) 



and p v 0 and p y correspond to the tem- 
peratures tg and t a d , respectively, for a 
wet surface. Because X varies with t_, the 
values of t g at each point must he computed 
by trial . 

4. In order to determine the heat transferred 
• through the blade metal of the segment, plot 
(tg^ d •* b a ) against s/c for the chordwise 

extent of the heated surface and obtain the 
average value by graphical integration or by 
Simpson's rule. The heat transferred through 
the blade metal Qtrans 18 then 

t 

^trans = ^g,d " ^s^av ^g A 8 ( 24: ) 


5. Determine the energy change due to friction in 
the segment F per pound of gas flow (from 
reference 12 ) 


1-2 


r,m 

2g 


JL f 

' Dh 


(25) 


where 


f = 0.0056 + | (Rej^ h ) 


-0.32 


(26) 


ana 


Re = w Ph 

^p hg^m (3600 g) 


(27) 


6. Determine the polytropic specific heat of 


gas c n from 


^ 1-2 = 


_ ^trans *1-2 


w 


- °n (*g,l - *8,8) < 28 > 
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and the polytropic exponent n from 


' c - c 
n = _2 JL 


(29) 


7. Determine the pressure of the gas leaving the 
segment p g ^g from 


•P, 


8,2 ~ P g, 


^(■v) 


n 

1-n 


(30) 


where ^ is estimated for a particular hot- 

gas-flow system from published experimental 
pressure-loss data on ducts and manifolds. 

Find, the change in radial kinetic energy E^.g 

per pound of gas flow through the segment from 


E- 


1-2 = 


u 2 - u 2 
r,2 r.l 

w 2 

2g 

~ (3600) 2 (2g3 


/ , x2l 

/ i Y - 

( I ^ 

V A P,2 P g,2/ 

VV 1 p s,v 


(31) 


8. Determine the net compression vork W 1-2 

(appendix D) done on each pound of gas in flow- 
ing through the segment from 


u t £ - u t y 

*1-2 - "prop " F 3>2 - h - 2 g ^ - 


" F l-2 


(32) 


9. Finally, check. the accuracy. of the assumed tg^g 
(appendix D) from 


: g,l " -8,2 ~ 


(«i-2 - %.2> <“) 


and reestimate . t g g until the value obtained 

from equation (33) agrees with the initial 
estimate in step III-A. 
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' B. For the subsequent radial segments, repeat step III -A 
using the values of t g ^ 2 and p g ^ 2 °f one segment 

as tg^i and p g ^, respectively, of the next 

segment. 

C. If t 8 at any point (found in step III -A -3 for the 
correct t g ^ 2 ) proves to he less than 32° F, it can 

he increased by increasing the hot -gas flow w and 
repeating step III. Usually t s will be a minimum 

at the stagnation point and at a segment near the 
hub . 


IV. From the final conditions of the gas leaving the last 
segment., determine the tip nozzle area A n required to 

maintain the gas flow (for unity discharge coefficient) 
from ' 



V. 


Calculate the following: The required input from the 

heat -source Qj, excluding losses in the induction sys- 
tem prior to the initial point, from 


% - » ° P (* 8 ,i - T o> - j 


(36) 


the total heat added to the gas flow Qy from 


( 37 ) 


the heat escape at the tip nozzle from. 
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On* w ' c p (T gr f -V - (3«) 

\ : ■ ' ' - 

and the heat transferred through the entire heated sur- 
face of the blade ZQfc rans from 


^ = * Cp (Tg'irT^f > + J (!*£ 


2 2 \ 
- u t,i ) 
2g V 


(39) 


Thus, the effectiveness! of the propeller blade as a heat 

, . . . ^trans . 

exchanger is given by , ■■ — 


NUMERICAL EXAMPLES , 
Flight and Icing Conditions 


The procedure for determining the design requirements for a 
gas -heated propeller has been .applied to a theoretical hollow steel 
propeller having blades with : ci’ntral radial ribs, as shown in 
figure 1. The propeller -blade -form characteristics are given in 
figure 4 and other pertinent data in table I. An assumption has 
been made that heated air enters the blade shanks from a collector 
ring through suitable orifices and passes radially outward through 
only the leading-edge cavity in the blades, as illustrated in 
figure 1. The blade was arbitrarily divided for convenience into 
four radial segments each 1 foot, long starting at the 18-inch 
station. The radial stations for which the surface -temperature 
calculations were made were therefore at 24, 36, 48, and 60 inches. 
(See fig. 4.) For more accurate results calculations nearer the 
hub may be necessary. A greater number of shorter radial segments 
may also be needed in regions where; either the blade airfoil sec- 
tion or the internal flow area changes rapidly with radius . 

For purposes of comparison, the calculations have been made 
for the two following operating conditions through the same icing 
condition: 


Condition 

A 

B 

Pressure altitude, ft 

.18,000 

18, 000 

Ambient -air temperature, °F 

0 

0 

True flight speed,, mph 

400 

525 

Propelley speed, rpm ' : ' 

. 1430 

: 800 

Li^uidirv^itgr content, g/cu m 

04 

0,4 
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The value of liquid -water content chosen was based on the 
recommendations shown in figure 2. 

. For convenience in determining the initial internal gas condi- 
tions at the blade shanks, an arbitrary assumption was made that 
the internal- static pressure in the blade shanks would be the 
ambient atmospheric pressure plus three -fourths of the f light, 
dynamic pressure. The assumed hot-gas flow was varied at an. — 
initial temperature of 500° F until surface temperatures of 32° F 
or more were obtained over the heated surface of the blade. 


Results 

The general results are summarized as follows : 


Condition 

A 

B 

Gas flow per blade, lb/hr, 

450 

750 

Initial gas temperature, °F 

500 

500 

Final gas temperature, °F 

349.6 

. 320.8 

Heat-source input required per blade, Btu/hr 

50, 550 

86,680 

Required blade-tip nozzle area, sq ft 

0.00466 

0.0132 


The results of the calculations are given in more detail in 
table II and figures 5 to 10. For the four radial stations and 
for both simulated flight conditions, the chordwise variation of 
.local velocity ratio V^/V-g is shown in figure 5; the external 

heat -transfer coefficient h a is shown in figure 6; the local 
water interception M and evaporation M^. rates are shown in . 
figure 7; the local heat -transfer datum tender ature t a ^ is 
shown in figure 8; -and the external -surf ace temperature t a is 
shown in figure 9. , The radial distributions of internal gas tem- 
perature tg, internal gas pressure p g , internal heat -transfer • 

coefficient hg, and the external-surface temperature at the 

stagnation point t 8 .are shown in figure 10 for conditions A 
and B. 

The calculations of the external heat -transfer coefficients 
were made using equations (1), (3), and (12). Attempts were made 
to locate accurately the external boundary-layer transition 
points using the methods of references 13 and. 14, and the theoret- 
ical, transition points were found to be .located aft of the heated 



18 


NACA TN No. 1494 


areas. In four cases the velocity distributions (fig. 5) have 
fairly distinct peaks near the leading edge and an arbitrary 
assumption was made that the surface roughness due to the imping- 
ing water would cause transition at these minimum-pressure points. 
This assumption applies to the camber face at the 48-inch station 
for condition A and to the thrust face at the 60-inch station and 
the camber face at the 36- and 48-inch stations for condition B. 
The resulting chordwise distributions of h a> shown in figure 6 

for these cases, are similar to the' type shown on the left side of 
figure 11 and discussed in appendix A, inasmuch as in each case, 
the transition point occurs very close to the leading edge. 

1 ’ . fr « . • . 

The surface temperatures ; are obtained from equation (22), in 
which the external heat-treu^eb^^s^:. arbitrarily assumed 
to extend aft of the blade -rib center line for a distance on each 
f ace equal" to half the maximum blade thickness . This assumption 
is to account approximately for the' heat transfer to the trailing 
half of the blade, which in reality diminishes all the way to the 
trailing edge. An unpartitioned blade would not require such an 
approximation. . V'. ' ' 

DISCUSSION 

Thin-Skin Approximation 

The largest source of error in this analysis is thought to be 
the thin-skin approximation used in the solution of external sur- 
face temperatures. This approximation makes use of two related 
assumptions . First, the value of Hg depends on the temperature 
differential (tg,d “ t a ), which the temperature gradient 
through the blade metal Is ^assumed to be zero. , This gradient and 
the error involved is usually small, except at the leading-edge 
region, for any propeller blade made of steel or metal with equal 
or higher thermal conductivity. Second, the heat balance employed 
in equation (22). is based on the assumption that the total heat 
transfer from the internal gas to .the external air can be con- . 
sidered to be distributed uniformly and to pass only normally 
through suitable internal and external heat -transfer areas; that 
is, the iblade metal is assumed to be infinitely thin and surface- • 
wise heat conduction is neglected. The ratio of the two heat- 
transfer areas is then used to determine the ratio of the local 
internal and external heat -transfer rates. The error involved in 
this assumption is thought -to be large inasmuch as the thermal 
conductivity of the blade metal is usually very high compared 
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with the conductivity through the gas and air surface films. Sur- 
facewise heat conduction within the "blade metal is certain to ; 
reduce the temperature gradients indicated "by the variations of 
surface-film heat-transfer coefficients and changes in metal 
thickness. 

The two foregoing errors can "be eliminated from the analysis 
if the blade segments are studied using the "relaxation" method 
developed in reference. 15. No attempt has been made to use this 
method because of the prodigious amount of time and effort involved. 
The method and a sample of the results to be expected from it are 
given in appendix E. ; 

Principal .Factors Affecting Heat Requirements 

From a consideration of the foregoing analysis and the numer- 
ical examples, the required internal heat. flow through hollow pro- 
pellers depends on the following factors exterml to the blade 
surface, which appear to have the largest effects in determining 
the local heat transfer through the metal necessary to prevent ice. 

Ambient temperature. - With other factors remaining constant, 
the heat transfer through the blade surface required for ice pre- 
vention as defined herein varies directly with the temperature 
difference between the freezing temperature and the datum tempera- 
ture (ambient temperature plus kinetic increment). Even though 
the assumed liquid-water content of the ambient air decreases with 
decreasing temperatures, the heat required to maintain a minimum 
surface temperature of 32° F steadily increases as ambient tempera- 
ture decreases; but, computed on the basis of immediate evapora- 
tion of all the water that strikes the surface, the heat required 
will decrease with decreasing temperatures because the liquid- 
water content decreases. The point at which the heat requirements 
are equal when computed on' both bases of heating represents the 
maximum possible value of heat required. This maximum value rep- 
resents excessive heating for the conditions used in the examples 
presented. .f 

. i" ' ' 

Water impingement. - An increase in the rate of water inter- 
ception requires a direct increase in the amount of heat required 
to raise its temperature above 32° F. A more important factor, 
when external heat requirements are considered, is that an increase 
in the rate of water interception increases the amount of runback, 
which increases the area of wetted surface and thus requires more 
heat to offset. the evaporative heat loss. The difference in tem- 
peratures between a - wet and a dry surface resulting from the same 
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internal heat flow can be observed in figure 9(b), where the 
surface becomes dry on the thrust face of both the 48- and 60-inch 
stations. The surface -temperature rise above ambient temperature 
more than doubles its value in the two. instances after the surface 
becomes dry. 

Kinetic heating. - Kinetic heating always reduces the heat 
requirements for ice prevention because it increases the heat- 
transfer datum temperature above the ambient temperature. Its 
temperature rise increases as the square of the resultant velocity 
and has a lower value in wet air than in dry air. This increment 
of temperature becomes important as a saving in required heat flow 
at high resultant velocities of, the- order, shown; in the two examples 
presented. A fact of interest, ; is; that heat requirements are 
reduced according to the square of the; resultant velocity (because 
of kinetic heating) and are also increased according to the square 
root of the resultant velocity, (because of the heat -transfer 
coefficient). The net result is; that as velocity is increased 
the heat required to maintain a Surface temperature of 32° F first 
increases and then decreases ; the location of the maximum heat 
requirement depends on the conditions involved, This maximum heat 
requirement probably occurs below a velocity of 400 feet per sec- 
ond for most conditions. Thus, except for other variations (such 
as changes in propeller -blade sections), the heat requirements for 
propellers on high-speed airplanes are usually greatest at the 
inboard stations . For the same reason, the heat requirements are 
reduced by increasing the rotational speed of the propeller. 

Location of transition. - The determination of heat transfer 
from a propeller blade depends to a large extent on the exact 
location of transition. As shown in figure 11, considerably more 
heat transfer results over the biade faces for which turbulent 
flow is assumed (as shown by the higher heat -transfer coefficients) 
than would be the case if laminar flow were assumed. The same 
result is shown in figure 9, in- which the low surface temperatures 
of the four faces under assumed turbulent flow indicate a higher 
rate of heat transfer. Of course, transition might be caused by 
water-film roughness on many of the surfaces for which laminar 
flow was assumed,' because the velocity gradients are quite low 
over all the stations, as shown in figure 5. f' : 

Altitude. - The altitude -pres sure term frequently appears in 
the analysis and its total effect on heat requirements is difficult 
to evaluate. The two principal effects of altitude (other factors 
remaining constant) on the required^blade-metal heat transfer are 
largely compensatory. With increasing altitude the heat transfer 
due to the evaporation processes increases whereas the heat transfer 
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due to external convection decreases'. Increasing altitudes have a 
serious effect on restricting the internal gas flow. In order to 
offset the reduction in gas density with altitude, the gas flow w 
must he maintained hy increasing the internal gas velocity ^ with 
resulting increases in flow friction and pressure drop. These 
increases, coupled with the loss of ram pressure at altitude, indi- 
cate that the required internal gas flow steadily becomes more 
difficult to maintain with increasing altitude and that a booster 
puinp may be necessary at high altitudes to attain the required gas 
flow. A pressure boost can also be. obtained by increasing the pro- 
peller, rotational speed .whenever possible. 


Comparison of Results of Numerical Examples 

. . Flight condition B required a larger internal heat flow than 

condition A for two reasons: the kinetic heating was less, owing 

to a lower airspeed and propeller speed, and the surface area 
under assumed turbulent flow was greater. ■ 

The required tip nozzle area for condition B was almost as 
large as the internal flow area at the exit to station 60, which, 
indicated that the gas., pressure at the tip had dropped nearly to 
its lowest limit and that a higher rate, of flow was almost 
impossible without augmenting the initial gas pressure or increas- 
ing the propeller speed. If this tip nozzle area is used, the 
flow will become excessive under conditions of higher airspeed 
and propeller speed or lower altitude. The hot.-gas flow could, 
however, be regulated by a throttle in the intake system. 

The required heat-source capacities for these two conditions 
appeared to be very high compared with external electrically 
heated blade shoes based on the meager experimental data avail- 
able. The high required heat-source input is, however, explained 
by the fact that internal gas -heating systems normally have no 
selective control of heat transfer and as a result much of the 
surface area is considerably overheated and further wastage 
occurs at the. tip nozzle. In the examples that used this method 
of blade heating, only £ narrow band at the leading-edge surface 
had temperatures near 32° F. Consequently, if a narrow strip of 
ice 0.3 inch wide, for example, were to be permitted along the 
leading edge, the internal heat -flow requirement for station 24, 
condition B, (most, critical) would be nearly halved. This 
reduction would result because ^he s^rfa^e temperature at the 
stagnation point could thep be allowed -to drop approximately . 

9° below 32° F and the differential |t a - t a a) would be 
approximately one-half its prigin^! va^ue ^17! 6° F). Because 
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the stagnation point is the controlling factor for the segment, 
the internal heat flow could then he reduced by nearly half and 
the blade-surface temperatures would remain above 32° F except for 
the narrow leading-edge band. In practice, the internal heat flow 
required therefore depends largely on whether narrow strips of ice 
along the leading edges are tolerated. 


Suggested Technique for Reducing Heat Flows 

l ' • • 

A method of reducing the. required heat flow in a typical hol- 
low propeller blade is suggested in which the internal flow passage 
is altered to attain a more economical- distribution of heat transfer. 
When it is assumed that equation i (18) t dhd' the hydraulic diameter 
concept can be applied to internal passages whose perimeters con- 
tain convolutions, designs >bf; fiow passages' can be developed whereby 
large reductions in heat-source capacity, and pumping-power output 
are apparent over the typical hollow-blade interior, as. shown in 
figure 1. Two suggested alterations are' shown in figure 12 and 
compared with the original blade section.. Using primes to denote 
the altered designs, the -following mathematical comparisons can 
be shown: 

' i 

For both the original and altered designs 1 

4.1 X 1CT 4 (T„ J 0 - 3 w 0 - 8 P 0 - 2 k e 

. ' ' Vg* 


from . equat ion (18). For the. same inlet gas temperatures, assuming 
P = Ag. and taking the ratio of a modified design to the original 

design 

? (Sl\ 

h g A g < w \A g ) \^<j 

Setting' this ratio equal to unity will permit evaluation of the 
flow ratio -that will produce, the same amount of blade-surface heat 
transfer as the priginal design, or , 
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** = ApPg(3600) 


ss 



Modifica- 
tion 1 

Modifica- 
tion 2 

0.149 

■ 

0.152 

! 

.298 , . 

.228 

I 

/ 

.813 

/ V . 

,481 


In s unmar i zing these modifications, if the area of the -inter- 
nal flow passage is reduced- and the internal heat-traftsfer area is 
increased, the same total heat transfer to the external surface 
can be maintained while the internal -gas flow is greatly reduced. 

A reduction of the flow lowers the internal velocity and the fric- 
tion loss and thus raises the available internal pressure. When 
the flow is reduced, more heat is" removed per pound of gas flow 
and, consequently, the gas temperature falls more rapidly and 
reduces the heat loss at the tip nozzle. Care should be taken 
with the internal heat -transfer area that some segments do not 
abstract such heat as- to starve subsequent segments. 

In addition to the foregoing reduction of heat flow, a better 
distribution of chordwise surface heat transfer can be accomplished 
by making the greatest increase of internal heat -transfer area at 
the leading edge .by the use of fins similar to those shown in 
figure 12. An accurate determination of the magnitude of the heat 
transfer resulting from these modifications requires that the 
whole problem be approached along the lines of the relaxation 
method as outlined in appendix E. All the indications shown lead 
to the belief that a hollow propeller blade can be designed with 
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modifications similar to those illustrated, which will efficiently 
provide ice protection, with a heat flow of feasible proportions. 

; COMMENTS AND RECOMMENDATIONS . 

The method presented for estimating heat requirements for ice 
prevention on gas-heated hollow, propeller blades provides inte- 
grated means for determining the design requirements for ice pre- 
vention on gas -heated hollow propeller blades operating in any 
selected condition. Several of the suggested simplifications and 
other short methods can be used where conditions do not require - a 
rigorous treatment. Solutions by trial are occasionally required, . 
which are admittedly inconvenient, but became fairly rapid with 
increased familiarity. i . 

Substantiation of several formulas , upon which this method 
is based, by experimental measurements in actual icing conditions 
is recommended. A relaxation analysis of the; whole; propeller 
blade should be made to determine the effect Of leading-edge 
internal ribs and the exact heat transfer normally, radially, and 
chordwise through the blade metal '* 

Flight Propulsion Research Laboratory, 

National Advisory Committee for Aeronautics, 

Cleveland, Ohio, June 27, 1947. 
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APPENDIX A 


EXTERNAL HEAT-TRANSFER COEFFICIENTS 

Airfoil coefficients based on boundary layer. - Several 
methods are available for calculating the local and average heat- 
transfer coefficients for the surfaces of airfoils (references 4, 
16, 17, and is). Local surface yaiues for the heat -transfer 
coefficients on airfoils (equation (l)) are derived in reference 4, 
primarily on the basis of Reynolds' analogy between skin friction 
and heat transfer through laminar boundary layers, and an alter- 
nate form of Reynolds ' analogy (equation- (3) ) is presented for the 
turbulent boundary-layer lease ^ in ref erenbe 16. These solutions 
are detailed because they- take, into--, account the variation of local 
air velocity and boundary-layer thickness over the airfoil surface. 
No account has been taken either in this' analysis or in the refer- 
ences, however, of the unknown effect on .the boundary-layer heat 
transfer of the heat addition ; to the boundary layer caused by a 
heated surface. - t ‘ : ' . •- • 

Empirical airfoil coefficients. - The methods of reference 8 
( equations (9) and (lO)) have been added as alternate solutions 
for h a because they offer easier solutions by their use of 
average air velocities over the faces of the airfoils. These 
equations require a slight but adequate approximation for use- in 
this analysis when determining the average temperature 7y, 

because the surface temperatures are unknown. However, estimates 
of average surface temperatures involving errors of less than 

1 « 

lg percent of h a can easily be made. 

Leading -edge cylinder coefficients. - At the stagnation point 
the value of h a cannot be determined in equations (l), (3), (9), 

and (10). The suggestion is made in reference 8 that the leading 
edge of an airfoil be regarded as an isolated cylinder in trans- 
verse flow for which heat-transfer data are available in refer-, 
ence 1C . When these data are used and Nusselt • s number at the 
stagnation point is expressed as a function of Prandtl number and 
Reynolds number, an empirical equation for the heat-transfer coef- 
ficient over the forward half of the equivalent leading-edge 
cylinder (equation (12)) is developed (reference 8). Figure 11 
diagrammatically illustrates the suggested technique of fairing 
the distribution of heat-transfer coefficients over the leading- 
edge cylinder into the laminar- and turbulent-flow distributions 
for the rest of the airfoil in two typical configurations. On tjip 
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thrust face, transition is shown to occur qlose to the leading 
edge and the h a curve is faired directly from the leading-edge 

values into the turbulent -flow curve. On the camber face, transi- 
tion is shown to’ occur after a length of laminar flow and the 
curve is arbitrarily faired into the turbulent regime. The former 
type of distribution occurs in the numerical example, in certain 
cases on both faces of the blade. 

Location of transition. ' - In the foregoing considerations, 
the location of the transition point must be determined. For 
increasing local pressures, the point of transition, which may be 
assumed to be coincident with lamina!* separation, maybe estimated 
by the method of references 13 and 14. For decreasing pressures, 
the suggestion is made in reference 5 that the Reynolds number 
based on the laminar boundary layer thickness Reg ^ at the point 

of transition is between the limits of 8000 and 9500. 

The location of the transition point greatly affects the heat 
transfer from an airfoil because the heat-transfer coefficient for 
turbulent flow is considerably larger than that for laminar flow. 
No reliable method is known for predicting the point of transition 
when the airfoil is undergoing water impingement. The presence of 
water on the surface probably causes transition to occur forward 
of the point determined by the foregoing methods. 
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APPENDIX B„. 

EFFECTS OF WETTED SUEFACE 

Rate of water interception. - The rate of water interception 
is only slightly overestimated by assuming that the droplets 
travel in straight-line paths intersecting the blade surfaces. 

This condition exists when the airfoil collection efficiency 
becomes 100 percent. Exact collection efficiencies for airfoils 
are unavailable, although fairly complete data are available for 
cylinders. For symmetrical and low-camber airfoils the use of 
existing data for cylinders whose diameters are equal or related 
to the leading-edge diameters of the airfoils can be made with a 
fair degree of approximation. Accordingly, prope 1 ler -b lade -a ec t i on 
leading-edge cylinders, under ordinary flight conditions, have col- 
lection efficiencies varying approximately from 80 percent near* 
the hub to 98 percent near the tip. The error in assuming : 

100-percent collection efficiency for the propeller sections is 
not so great in figuring the rate of water interception as .in 
determining the extent of jetted surfaces, as discussed in the 
' succeeding section. 

The local rate of water interception M is thus proportional 
to the sine of the angle ^ subtended by the relative wind 
(helical speed) vector and the tangent drawn to the blade surface 
at the point being considered. The value of M is given by 
equation (13). 

Surface evaporation. - The effect of evaporation of water 
from the surface of -an airfoil in increasing the heat transfer - 
was first derived by Hardy (reference 19). He has presented an 
expression for a local evaporation factor X, which multiplies 
the local convective heat-transfer coefficient wherever the sur- 
face is wetted. This factor is given by equation (23) and the 
local rate of evaporation of water from the surface is given by 

'Hot - ^ (X-D (t a - t a ,d) <4Q\ 

Xj 

Where the local rate of evaporation is larger than the loc$!\. 
rate of water interception, the surface will tend to become dry, 

If the blow-off of water from the surface is neglected, however, 
there will "be runback of unevaporated water from the leading-edge 
region when the rate of intercepted water exceeds the rate qf 
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evaporation. In this analysis the surface of the propeller-blade 
section is considered dry .aft of the point s/c at which 

’'s/c I ps/c 

M d(s/c) ‘= ! M ev d(s/c) 

>-0 , ■ v . 0 

where the lower limit 0: is at the stagnation- point. For a dry 
surface, X becomes 1 and "no surface evaporation occurs. 

: ; , •* . ■ 

Heat -transfer datum; air temperatures - The heat -transfer 
datum temperature t a at j any : point, on ; the. blade surface is the 

■} ■ ,■ J’ J 11:' , f < ; i*' i . i- 

surface temperature that • would be obtained if . the blade were a 
nonconducting unheated bbdy. The; yalub \of ; ,;t a d is given by 

equation (14) based on figure. 3 and i equations; (16) and (17). The 
nature of t a d for a typical propeller; blade operating in either 

saturated or unsaturated air is shown in ; figure 13. There -are two 
components of t a ^ d : the adiabatic temperature t t or t t ^ w at 

the outer edge of the boundary layer and the temperature differ- 
ential At a or At a w in the boundary layer due to kinetic fric- 
tion. The frictional temperature rise in the boundary layer is. 
slightly larger in turbulent flow than in laminar flow according 
to equations (16) and (17) (as the Prandtl number is raised to the 
one-third power in turbulent flow compared with the one-half power 
for laminar flow) . Both the components of t & d are affected by 

the presence of water. The frictional temperature rise in a 
boundary layer that remains saturated is less than the rise in an 
unsaturated boundary layer, due to evaporation. Likewise, the tem- 
perature at the edge of the boundary layer for a saturated air 
stream follows the saturated air adiabatic line. 

Whether or not the air is saturated at a given point is often 
difficult to determine and, consequently, the determination of 
t a d is uncertain. In this analysis the following course has 

been adopted: For a wet surface, which exists forward of the point 

where the accumulated water intercepted equals the accumulated 
water evaporated (as defined in the previous section) 

ta,d = ^b,w + A1;i a,w 

For a dry surface, which is assumed to begin abruptly aft of this 
point 
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V = V + At a (4l) 

Accordingly, only when the ambient air is unsaturated is 

*a,d = + At a 
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APPENDIX C 


HEAT BALANCE 

Local Internal heat transfer. - The internal gas -film heat- 
transfer coefficient h_ is given by equation (18), which appears 

O 

to be the most applicable fonnula for determining this coefficient. 
The equation is based on data obtained for fully developed turbu- 
lent flow in long tubes. It is doubtful, however, whether this 
equation will give satisfactorily accurate results when applied to 
flow through hollow propellers for two reasons: (1) The nonsym- 

metrical cross sections of a hollow propeller-blade passage and 
the effect on heat transfer of changes in flow area may not be , 
fully accounted for in an equation based on flow in straight cir- 
cular tubes; and (2) no account is taken in equation (18) of local 
flow variations along a plane normal to the direction of flow 
(chordwise in a propeller blade). In a rotating propeller, the 
Coriolis acceleration produces radial velocities along the leading 
edge of the internal passage that are higher than those toward the 
trailing edge. The local internal heat-transfer coefficient would 
therefore probably be higher at the leading edge than is indicated 
by the average value given by equation (18). 

When the value of h g given by equation (18) is used, the 
local internal heat transfer for a unit area can 5 be written as 

°6 * h 6 <* 6,4 - *-> (43) 

\ 

where t g ^ is the mean-flow static gas temperature plus the tem- 
perature rise due to surface friction (equation (20)). 

Local external heat transfer. - The local rate of external 
heat transfer per unit area at any point may be expressed as 


H a = h a X (t 8 - t a ^ d ) + M Cjh (t 8 t Q ) - ?LgB_ (44) 

where c m is the heat capacity of liquid water, taken herein as 
unity. The first term of equation (44) on the right-hand side 
gives the combined convective and evaporative heat transfer; the 
second term gives the sensible heat absorbed by the intercepted 
water during its temperature rise to the surface temperature; and 
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the third tern gives the heat released by the impinging water at 
the expense of its kinetic energy. For simplicity, an approximate 
expression for H a can be used, wherein 

I a = (h a X + M) (t s - t a<d ) (45) 

and the error in the approximation is small except at the leading 
edge where M is a maximum. 

Heat-balance- equation. -In this analysis the equilibrium of 
heat transfer. to and from the. blade surface at any point is 
expressed as . 

. ' H g = H a (46) 

^a 


When equations (43) and (44) are substituted in equation (46) 
t g can be solved for: 


t D = 


^6^8, d A g + 


I hoXt, 


+ M ( t o . + §j). 


Aa 


hgAg + (h a X + M) A fl 


( 22 ) 


Equation (46) is based on the assumption that a thin-skin 
approximation can be used for determining the surface temperatures. 
(See section entitled "DISCUSSION.") 
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. APPENDIX D 

THERMODYNAMIC ANALYSIS OF GAS FLOW 

. - - ( t ■ ( 

Assumptions. - The method for determining the temperature and 
pressure of the internal gas during its flow through the hollow 
propeller blade is based on the thermodynamic steady -flow energy 
equation* The internal flow is assumed to have a uniform radial 
velocity distribution at any station. The assumption is made that 
for a short radial segment the polytropic exponent n of the 
internal flow process can be considered constant. The error 
involved will v^ry as the radial; '/length l of the segments used 
in the solution. For convenience in; using the polytropic rela- 
tions for gas, the energy due to flow friction is assumed to be 
heat added to the gas at the. expense of work done by the gas. 

Analysis of flow process. - The ' thermodynamic steady -flow 
equation may be written in terms of the energies per pound of gas 
as 


U “ -* u 

*1-2 - J Ql-2 = r ' 2 "2^ - " £ ' — + J ( U 2-V + (P2 V 2 ‘ Pl v l> ( 4T > 

This equation states that the difference between the net compres- 
sion work done on the gas and the net heat lost from the gas in 
passing t hr ough the segjnent is equal to the total of the differ- 
ences between the kinetic, internal, and potential energies at the 
leaving and entering ends of the radial segnent. 

For a mechanically reversible process, the sum of all the 
mechanical -energy terms in equation (47) can be expressed by 

E l-2 + (p 2 v 2 - Pl v l> " W l-2 < 48 > 

where E^.g is defined as the radial kinetic energy term by 

equation (31) , The net work for a polytropic process satisfying 
the relation* pv 11 = constant is given by 

<‘ 6,1 - (49) 


and the flow work can be expressed as 
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Vs - Vi = E(t e,2 - W < 50 > 

When equations (48) and (49) are equated and equation (50) is sub- 
stituted, the net compression work hecames 

w l-2*'®l-2.*^f»<*g,2 - ‘ 8 ,l) <51) 

\ - , 

In a similar manner, the internalrenergy term of equation (47) can 
be written as 

J < B 2 - U l > = £ <‘ 8,2 * (52) 

When equations (50), (51),. and (52) are substituted, equation (47) 
becomes 

\ 

Q l-2 = (lrn)(^l) J ^g,2 " ( 53 ) 

However, 

o - ^trans E l-2 ^ s 

Q l-2 ~ ~ - — = c n(%l - tg; 2 ) . (28) 

The gas -temperature differential through the segment can now be 
expressed in terms of the energy quantities by combining 
equations (51), (53), and (28) to obtain 


" *g,2 


7-n 

(y-l) uJc. 


n 


( W l-2 " E l-2^ 


(33) 


In equation (33) the value of W^_g is given by equation (32), 

The increment of enthalpy W pr0 p due to work done by the rotating 

propeller on a pound of gas in flowing through a radial segment is 
given by 


W 


prop 


^,2 


- 


2g 


(54) 


This equation can be shown as follows: The propeller-drive-shaft 

torque increment due to the internal gas flow through a blade 
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segment is given by the change in moment of tangential momentum, 
which may be written in terms of torque per pound of gas per 
second as 


* = *2^,2 ~ Tl^t.l ' - 

, 8 ' 


(55) 


where Uj. is the blade tangential._yelocity at the radius r. The 
rate of work done on a. pound of gas can be written as 



< ar 2 u t,2 


«t ) 2 2 - 

& , ' " g 



2 

1 


(56) 


because u t = tor. For a straight rotating* channel, the increase in 
tangential kinetic energy between two points is obviously given by 




u. 


,2 

hi 


for a pound of gas, so the remainder of the total work represents 
the amount of external work done on a pound of gas in flowing 
through a radial segment to produce a change in the enthalpy of the 
gas. 
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APPENDIX E 


EFFECT OF HEAT CONDUCTION THROUGH BLADE 
METAL AT LEADING EDGE 

A study was made using the relaxation method as described in 
reference 16 to determine the effect of heat conduction in the 
blade. skin on the temperature distribution around the leading edge. 
The original assumption that heat is transferred only in a direc- 
tion normal to the blade skin, which was assumed to be infinitely 
thin, was expected po be appreciably upset by the variations of 
blade-metal thickness and of the effective external heat transfer 
around the leading edge. 


Laplace ' 8 equation, which requires the reasonable assumption 
that the blade metal is a homogeneous, isotropic solid, was used 
with conditions at the internal boundary defined by 


k = h i 
dy S 


(57) 


and at the external boundary by 

d0 

k 35 : = (h a X + M)(3 - 0) 

where 


k 


0 


t 


X 


y 


g 


(h a x + m) 


thermal conductivity of blade metal, 
Btu/(hr)(sq ft)(°F/ft) 

<*g,d - *,). *> • 

heat -transfer datum gas temperature, °F 

metal temperature (surface or. internal), °F. 

distance normal to boundary ft 

internal gas heat-transf er coefficient, 
Btu/(hr)(sq ft)(°F) 

effective external heat -transfer coefficient, 
Btu/(hr)(sq ft)(°F) 


(58) 
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<* 8 ,d 

o 

heat-transfer datum air temperature, F 

" The scale used vas 1 inch = 0.05 inch and for convenience the 
solution was arbitrarily originated 1 inch in chord length from 
the leading edge. The metal temperatures at. the starting points 
on both faces of . the blade were, assumed to be slightly less than 
the values of surface temperatures that were obtained in the numer- 
ical example. The inclinations ( af the isotherms to the surface at 
the starting points were estimated from the anticipated nature of 
the heat transfer. The values of metal temperature obtained are 
unreliable, but the isotherm and surface -temperature patterns 
obtained are believed indicative of the trends to be expected if 
the solution were to be extended over the entire blade cross 
section. : .■■■■. '■ 

The calculated results ^obtaihed for. simulated flight and 
icing condition B at the 24-inch’station are shown in figure 14. 
Because heat is transferred everywhere within the metal normal 
to the isotherms, a considerable' quantity of heat flows through 
the metal of the camber and thrust faces into the mass of metal 
at the leading edge, as shown in figure 14. The surface tem- 
perature distribution is changed and thus the surfacewise tem- 
perature gradients at the leading edge are greatly reduced. The 
surface temperature rise at the leading edge due to conduction in 
the metal is approximately as indicated by a comparison of the 
two curves at the bottom of figure 14. 
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TABLE I - PROPELLER DATA 


Station 

24 

36 

48 

60 

• 

Blade radius, r, ft 

2 

3 

4 

5 

Blade -segment radial 
length, l, ft 

1 

1 

1 

1 

Blade chord, c, ft 

1.05 

1.06 

1.07 

0.986 

Equivalent leading-edge 
cylinder diameter, 

'D c , ft 

0.027 

0.0095 

0.0057 

0.0040 

External heat -transfer 
area, A^, sq ft 

1.26 

1.18 

1.15 

1.06 

Internal heat -transfer • 
area, A g , sq ft . 

1.13 

1.05 

1.01 

0.942 

Internal "blade-passage 
perimeter, P, ft 

1.13 

1.05 

1.01 

0.942 

Internal blade-passage 
cross -sectional area, 
Ap, sq ft 

0.0558 

0.0315 

0.0239 

0.0173 

Internal blade passage 
cross-sectional area at ' 
point 1, Ap^, sq ft 

. 

0.0875 

0.0375 

0.0275 

0.0205 

Internal blade passage 
cross-sectional area at 
point 2, Ap 2 * Bq, ft 

0.0375 

0.0275 

0.0205 

0.0141 
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Figure 3. - Adiabatic temperature and pressure chart for saturated and unsat- 
urated air- (Data from reference 9.) 




NATIONAL AOVISORY 
COMMITTEE FOR AERONAUTICS 


g. 4 


NACA TN NO 













Velocity patio, V^/V 


NACA TN No. 1494 


Fig. 5 





(a) Condition A. 



(b) Condition B. 

• - Chordwise variation of local velocity ratio for two con- 

d i t ions. 


Figure 5 











External heat-transfer coefficient, ha, Btu/(hr)(sq ft)(°P) 


Fig. 6 a, b 
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Fig. 6 c, d 
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Fig. 7 a, b 
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Figure 7. - Chordwise variation of rates of water interception and water 

evaporation for two conditions. 
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Fig. 




Figure 7. - Concluded. Chordwise variation of rates of water inte 
ception and water evaporation for two conditions. 
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Figure 8 



External-surface temperature 



Figure 9 
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(a) Condition A. 

Figure 10. - Radial variation of internal gas temperature, internal gas press- 
ure, internal heat-transfer coefficient, and external-surface temperature at 
stagnation point. 
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Figure 10. Concluded* Radial variation of internal gas temperature, 
ternal gas pressure, internal heat-t ransfer coefficient, and external 
surface temperature at stagnation point. 




Diagram of typical chordwise variation of external heat-t ransfe r 
coefficient over propeller blade. 
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Fig. 14 



Figure 14. - Comparison of leading-edge surface temperature by two meth- 
ods for station 24, condition B. 







